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Abstract. The earth tide provides a spatially 
and temporally predictable force that deforms the 
Earth and can be measured as changes in gravity, 
tilt, and strain at the Earth's surface. All 
things being equal, tidal constituent amplitudes 
and phases will not change with time. However, in 
the vicinity of earthquake focal regions con- 
ditions may not be equal. Crustal rocks stressed 
to more than ~0.6 of their failure strength ex- 
hibit material properties over and above that of 
linear elasticity. These effects, including 
dilatancy, are known from laboratory measurements 
but have not been proven situ . 

interactions between the earth tide and crustal 
rocks that are under high tectonic stress are 
discussed in terms of simple phenomenological 
models. In particular, the difference between a 
nonlinear elastic model of dilatancy and a 
dilatancy model that exhibits hysteresis is noted. 
It is concluded that the small changes in stress 
produced by the earth tide act as a 'probe' of 
the properties of crustal rocks . Observations of 
earth tide tilts and strains in such high stress 
zones may, therefore, provide keys to the con- 
stitutive properties and the tectonic stress rate 
tensor of these zones. 

Introduction 


Research concerning earth tides has concen- 
trated almost exclusively on the elastic behaviour 
of the Earth's crust and mantle, and the bound- 
aries between elastic materials that give rise to 
the so called geologic, topographic and cavity 
effects. There are, however, zones of high 
tectonic stress within the Earth's crust that may 
have a more complex rheology which can be invest- 
igated with the earth tide. Small changes in 
stress associated with the tide act as a 'probe' 
both of the state of tectonic stress and of the 
constitutive properties of these zones. This 
concept, which has yet to be observed in iii situ 
experiments, is based on results from laboratory 
measur^ents of the continuum properties of 
intact rock samples subject to deviatoric stress 
in excess of ~0.6 of their failure strength. 

These ideas are discussed in terms of simple phen- 
omenological models. A rigorous treatment, 
including other possible rheologies, is not 
warranted until observations exhibiting anomalous 
tidal variations have been recorded. 

Beaumont and Berger (1974) rather naively 
suggested that temporal variations in the Earth’s 
admittance to the earth tide should accompany 


V /V seismic velocity anomalies, if the seismic 
p s 


velocity anomalies result from changes in the 


elastic properties of the crust. Their theoret- 


ical predictions suggested significant tidal 
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admittance anomalies in the vicinity of 'elastic 
dilatant' crustal inclusions. These results have 
been criticized from the standpoint that the 
dilatant behaviour of rocks as observed in the 
laboratory {see, for example, Brace et al, 1966) 
is not merely a change in linear elastic proper- 
ties. Consequently, an elastic inclusion model 
for tidal response represents an oversimplific- 
ation. This paper provides an opportunity to 
extend Beaumont and Berger ' s results to more 
realistic rheologies. 

Tidal Interactions with a Nonlinearly Elastic 
Crust 

In this paper the normal definition of dilatancy, 
a volumetric change induced by a deviatoric stress, 
is used. Infinitesimal linear elasticity theory 
is represented by the tensor equation 
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where £ . . and a. . are the strain and stress 
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tensors, 6^^ is the Kronecker delta, A and G are 

constants and Ij is the first invariant of the 
stress tensor = 0^^ (surimiation convention over 

repeated indices implied) . This equation is 
separable into isotropic and deviatoric parts. 
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which demonstrates that deviatoric stress, 

D 1 

O . . = O. . - -rO. . , can induce only deviatoric 
13 13 3 11 ^ 

strain and no volumetric change. Therefore, 
linearly elastic materials undergoing infinites- 
imal strain do not exhibit dilatancy. 

Stuart and Dieterich (1974) considered a model 
of nonlinear elasticity complete to second order 
in stress (Reiner, 1945) , 

= <(>o 6 ^^ + + (j> 20 ^^ , i,j = 1,2,3 , 

where cj>j^ = »iz < I 3 ) k = 0 , 1,2 , 

and Ii , 1 2 and 1 3 are the scalar invariants of 
the stress tensor. Such a model requires seven 
instead of two elastic constants, that is a, B, 

C, H and M in addition to A and G. 

e. . = (a+Aii+Bii+ci 2 ) 6 . . + (G-Hiii)CT. ,+Ma., a, . , 
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i,j,k = 1,2,3 . 

These constants determine the quadratic 
behaviour of nonlinear elasticity. Isotropic and 
deviatoric stress-strain relationships are cross- 
coupled, therefore, deviatoric stress induces 
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volumetric changes. In addition, the behaviour of 
such nonlinear elastic materials is anisotropic. 
This anisotropy, as emphasized by Stuart and 
Dieterich, is not intrinsic. It is a physical 
anisotropy that is caused by the stresses, 'but 
not the stress-strain relationship itself! 

A possible nonlinear stress-strain relation- 
ship (figure 1) illustrates the interaction among 
a nonlinear elastic material, the tectonic stress, 
a, and the earth tide stress, a^. This simplified 

approach ignores the fact that the a versus e or 

0 curves are functions of the stress invar icints. 

No matter what the state of tectonic stress, 
strain in any direction will follow a curve of 
the type shown in figure 1, either concave or 
convex. The tidal stress (~0.01 bar) , being much 
smaller than the range of tectonic stress (~5 
kbar) or earthquake stress drops (~100 bar) , will 
sample an essentially linear elastic behaviour 
proportional to the elastic teingent moduli of 
the nonlinear tectonic stress-strain curve. 
Therefore, the tidal admittance remains linearly 
elastic, but almost certainly anisotropic, and 
changes with increasing tectonic stress. Figure 

1 suggests a decrease in tidal strain with 
increasing tectonic stress, but the exact form of 
the change will be more complex, some components 
of the strain tensor may increase while others 
decrease. Such changes can in theory be pre- 
dicted in detail from the constitutive law and a 
knowledge of the tectonic and tidal stress ten- 

Model of Dilotoncy 



Fig. 1. A typical stress (a) versus strain (e) 
or dilatation (0) graph for nonlinear ly elastic 
crustal rocks. Superimposed tidal stress (O^) 

induces a tidal strain (£^) proportional to the 

tangent to the tectonic stress-strain curve. 


sors. At present,- mere identification of temporal 
variations in the tidal admittance would be 
sufficient to demonstrate both a rheology more 
complex than linear elasticity and a changing 
tectonic stress. Heaton (1975) has outlined some 
of these ideas and has also pointed out that anis- 
otropy will in general produce phase changes in 
the tidal admittance in addition to the amplitude 
variations. 

Tidal Interactions with a Crust that 
Exhibits Stress Hysteresis 

Intact rock sampies* when subject to cyclical 
stress that approaches the failure strength of the 
rock, exhibit dilatancy in the manner predicted by 
the Stuart and Dieterich model. That is, the 
stress-strain relationship is no longer linear for 
stresses in excess of ~0.6 of the failure strength 
and deviator ic stress induces volumetric strain. 
However, such samples exhibit stress hysteresis in 
addition to a nonlinear behaviour. Strain is a 
double valued function of stress, the value 
depending on v;hether stress is increasing or 
decreasing. Results from Scholz and Kranz (1974) 
(figure 2) on intact Westerly granite illustrate 
this behaviour. Repeated stress cycles (1, 18, 

19) reduce the minimum deviatoric stress required 
for the onset of dilatancy to ~300 bars, but in 
each stress cycle a broadly similar hysteresis 
loop is fillowed. A widely accepted explanation 
is that the energy loss, measured by the area of 
the hysteresis loop, is the work done against 
friction in opening and closing pre-existing 
microcracks within the rock sample (figure 2) . 

An idealized representation of the behaviour of , 
the 0 (dilatatipnal) , <}i and z components of strain 
in the cylindrical rock sample (figure 2, lower 


Dilotoncy of Intoct Westerly Gronite (Scholz ond Kronz, 1974) 



Fig. 2. Dilatancy of intact Westerly granite 
(from Scholz and Kranz, 1974) . The panels (left 
to right, top to bottom) illustrate: 1) the 
opening of a microcrack under applied stress; 2) 
and 3) the changes in dilatation (0) during 
stress cycles 1, 18 and 19; 4) idealized stress- 
strain relationships for the E^, 0, and E^ 

components of strain. 



Fig. 3. The variation of superimposed tidal and 
tectonic stress (O) with time (t) . “niis 
simplified graph illustrates a linearly 
increasing tectonic stress and a tidal stress 
comprised of a single sinusoid. The signific- 
ance of A, B and C is discussed in the text. 


panel) demonstrates that all strain components 
are double valued. This model assumes that all 
cracks start to open and , close at the same state 
of stress. It provides a first approximation 
to the true behaviour of the rock sample and can 
be represented mathematically in terms of an 
elastic-plastic material with internal back ■ 
stresses. This description is discussed in 
detail in a later section. 

The main concern is the interaction between the 
earth tide and crustal zones that exhibit such a 
dilatant behaviour with hysteresis. It will be 
seen that the interaction is different from that 
predicted for a nonlinearly elastic zone. Con- 
sider three simple cases (figure 3) in which the 
total stress is represented by a superposition of 
a linearly changing tectonic stress and a simple 
sinusoidally varying tidal stress. In A, the 
tectonic stress rate is either zero or very slowly 

varying (do /dt « (dCT /dt) ) . In B, the 
T £j m6dn 

tectonic stress rate is approximately equal to the 


mean tidal stress rate (dcr^/dt 


(da /dt) ) 

E mean 


In C, the tectonic stress rate is sufficiently high 
that there is no decrease in total stress (da /dt 
> (da 7dt) ) . 

E max 

The interaction of these three possible stress- 
time regimes with the idealized hysteresis curve 
(figure 4) suggests that the character of the tidal 
response (strain or tilt) is peculiar to each of 
the three tectonic stress rate regimes. Conse- 
quently, observations of tidal strain can, in 
principle, be used to measure tectonic stress 
rates . , 

In figure 4A the tectonic stress rate is much 
less than the mean tidal stress rate. For states 
of total stress , a, that are below that required 
for dilatancy (line 1-2) a superimposed tidal 
stress (a^) produces a normal amplitude elastic 



Fig. 4. Interaction between the earth tide (a^, 

e„) and an idealized stress hysteresis loop for 
E 

the three tectonic stress rates. A, B, and C 
from figure 3. 

tidal strain (e ) . Moreover, the same response 
E 

occurs even when the total stress is sufficiently 
high to produce a dilatant response (line 2-3) . 
This result is a direct consequence of the general 
character of hysteresis curves. That is, any 
small stress decrease will be accompanied by an 
elastic strain recovery. The dilatant strain is 
not recovered .until all the elastic strain has 
been recovered. Consequently,- for small tectonic 
stress rates the crust retains its normal elastic 
admittance to the earth tide stress . 

In figure 4C the tectonic stress rate is 
sufficiently high that there is never a decrease 
in total stress. In the pre-dilatancy region the 
response is identical to that of A. There will 
be no tidal anomaly. However, once the tectonic 
stress has moved into the dilatant region the 
strain is forced to change, in proportion to the 
slope of the 2-3 line because there is no stress 
recovery. The tidal strain undergoes an anomalous 
decrease, or increase (see figure 2) but remains 
linearly related to the tidal potential. 

The most interesting interaction is that shown 
in figure 4B. In the dilatant region stress 
reduction is accompanied by an elastic strain 
recovery proportional to the slope of line 1-2. 

The response remains elastic until the stress 
increases to the former stress maximum. Stress 
increases beyond the previous maximum induce 
strains proportional to the slope of line 2-3. 

The overall response (strain or tilt) is non- 
linear. The magnitude of the nonlinearity de- 
pends on the relative slopes of lines 1-2 and 
2-3, and on the relative size of the tectonic and 
mean tidal stress rates, the maximum nonlinearity 
occurring when these are equal. The main conse- 
quence of a nonlinear response will be the 
appearance of numerous additional lines in the 
tidal spectrum at sum and difference frequencies 
of the tidal constituents. Their detectability 
depends on the signal to noise ratio at the non- 
linear interaction frequencies. 
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Fig. 5. Elastic-plastic model of a weak crustal 
inclusion. The model has uniform elastic prop- 
erties (A and U) but the yield strength, 
characterized by ki , ka,. and ka, varies with 
position. Details of the precise yield criteria 
are discussed in the text. 

Plastic Models of Crustal Inclusions that 

Exhibit Hysteresis 

It is anticipated that only confined volumes 
(crustal inclusions) will experience high stress 
concentrations at any given time. Only these 
inclusions will have a rheology that deviates from 
the normal elastic response. Alternatively, 
inclusions that have repeatedly experienced fail- 
ure will almost certianly appear 'weak' in the 
sense that they will exhibit a departure from 
linear elasticity at lower stresses than the 
surrounding intact crust. An inclusion model 
(figure 5) similar to that employed by Beaumont 
and Berger (1974) but with an elastic-plastic 
rheology is used to explore the latter alternative. 
The model, a part of an cixisymmetric half space 
with a 10 )on radius, 10 km deep disc shaped 
inclusion buried at a depth of 4 km, is intended 
to be very simple. The half space has uniform 
elastic properties (A ,p ) but the inclusion is 
characterized by yield strengths ks and ka that 
are less than that, ki , for the surrounding crust. 
The result of particular interest is the 
character of the surficial strain, e(r), or tilt 
as a function of tectonic stress, O^, The results 

(figures 6 and 7) were obtained for two different 
yield criteria and associated plastic flow laws 
by finite element modelling (Bathe, Wilson and 
Iding, 1974, and Bathe, Ozdomir and Wilson, 1974). 

The radial strain anomaly as a function of 
increasing stress (figure 6) is for a model with 
a Von Mises failure criterion (Prager and Hodge, 
1951) . 

F(k,a) = ji _ k , 

where Jx is the first invariant of the deviatoric 
stress tensor and k = <7^/3, where is the yield 

stress in simple tension. Yielding occurs when 
F >. 0. The curves 1, 2 and 3 illustrate the 


Fig. 6. Surface radial strain anomaly (F(r)) as 
a function of three equal increases (1, 2, 3) of 
tectonic stress (^.j,) for the elastic-plastic 

model with inclusions that have a Von Mises 
yield criterion. See text for details. 

strain anomaly for three equal stress increments. 
The strain as a function of position is constant 
for stress states F < 0. However, once yielding 
has occurred (curve 1), the buried inclusion 
appears 'weaker' and a strain anomaly develops in 
response to increasing stress by virtue of the 
contrast in properties between the inclusion and 
the surrounding crust. The inset figures illus- 
trate the time variation of surface radial strain 
at selected points under the assumption that 
stress increases linearly with time. These 
figures may be interpreted in a similar manner to 
figure 4 because tectonic stress is proportional 
to time. The transition from an elastic (dashed 
line) to an anomalous response is now smooth, 
unlike that of figure 4 , because the surrounding 
elastic crust 'filters' the plastic response of 
the inclusion. Strain recovery (not shown) on 
stress reduction will follow a hysteresis curve 
very much like that observed by Scholz and Kranz 
(1974) (figure 2, smooth curves) . 

In fact, there is an almost perfect analogy 
between the behaviour of the elastic-plastic in- 
clusion model and that of a rock sample in the 
laboratory. For the rock sample, yielding (slip 
on pre-existing micro-cracks) occurs over a 
range of stress because each micro-crack has its 
cwn yield strength. Similarly, the crust is 
inhomogeneous and includes many inclusions like 
that of figure 5, each of which has its own 
yield strength. The elastic-plastic inclusion 
model is equivalent to a rock sample with one or 
more micro-cracks concentrated in one region. 

The analogy between a rock sample and the crust 
may be made conplete by choosing a plastic yield 
criterion which has the same form as that for 
slip across micro-cracks . A suitable criterion 
is the Mohr-Coulomb criterion, which when 
generalized to three dimensions becomes the 
Drucker-Prager yield criterion, 

F(Ii,J2,k>) = all + - k' 
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Fig. 7. Surface radial strain anomaly (e(r)) as 
a function of three equal increases (1, 2, 3) of 
tectonic stress (O^) for the elastic-plastic 

model with inclusions that have a Drucker-Prager 
(generalized Mohr-Coulomb) yield criterion. See 
text for details . 

(Drucker and Prager, 1952) , where Ja is the 
second invarieint of the deviatoric stress tensor 
and a and k' eire material properties related to 
the cohesion c and angle of friction l(J. 
a = 2sin\p//3(3-siniJj) , k' = 6ccosi{///3(3-sini|j) . 

It is also interesting that the associated flow 
law is dilatant. 

0^ = 3oX , 

«p 

where 0 is the rate of plastic dilatation and X 
is a constant proportional to the rate of plastic 

working, a. .e?., or defined in terms of the 
13 13 

plastic strain rate and the yield function. 

e?. = X3F/aa. . 

13 13 

(Drucker and Prager, 195 2) . 

The results for the Drucker-Prager model 
(figure 7) suggest that the surface tectonic 
strain anomalies will not be as large as those 
for the Von Mises model. The difference is partly 
a result of the two yield crieteria, though values 
for the constants a, k' and k were chosen so that 
yielding occurs at the same tectonic stress for 
both models. A more important factor is the 
effect of dilatancy. The major difference between 
the strain for the two models at distances ^IlO km 
is due to strain induced by tectonic uplift of 
the zone over the dilatant inclusion. As the 
value of a is reduced toward zero, the results of 
the Drucker-Prager model trend smoothly to those 
of the Von Mises model. 

When the tectonic stress is reduced, the strain 
behaviour of the inclusion models is very similar 
to that of rock samples. The plastic inclusions 
do not possess internal 'back stresses'; there- 
fore, the plastic strain would be irrecoverable 
if the inclusions were not embedded in an elastic 
matrix. 'Back stresses' in the elastic matrix 
ensure that the plastic strain is recovered. That 
there is irrecoverable plastic work done during 
this process ensiures that the stress-strain 


relation will exhibit hysteresis over a stress 
cycle. The situation in a stressed rock sample is 
very similar. The micro-cracks have no intrinsic 
'back stresses.' It is the stress in the elastic 
matrix that ensures that dilatancy is recovered. 

The only weakness of the plastic flew laws 
that have been used is that the rheology is 
elastic-perfectly plastic; that is, there is no 
strain hardening. Such a model suggests that 
once slip across a micro-crack has been initiated 
it will proceed at constant stress. 

We are now in a position to predict the form of 
tidal interactions with the elastic-plastic 
inclusion models. The e versus t inset graphs of 
figures 6 and 7 are the equivalents of the £ 
versus a graphs of figure 4. The interpretation 
is exactly the same with the addition that the 
tidal anomalies will vary with position on the 
surface of the model. At large distances from 
the inclusion the response remains normally 
elastic: no anomalous tidal admittance is pre- 
dicted. In the neighbourhood of the inclusion 
the character of the tidal admittance will 
depend on the relative tidal and tectonic stress 
rates. 

Discussion and Conclusions 

The intent of this paper has been a discussion 
of phenomenological models of linear and non- 
linear variations in tidal admittance that result 
from stress dependent properties of crustal rocks. 
The proposed constitutive relations are those 
observed for intact laboratory samples at devia- 
toric stress levels in excess of -0.6 of their 
failure strength. If crustal rocks exhibit the 
same properties situ , it follows that non- 
linear variations in tidal admittance of the kind 
predicted would indicate: 1) that tectonic 
stresses are sufficiently large to induce non- 
linear behaviour, 2) that crustal dilatancy with 
hysteresis is occurring, and 3) that tectonic 
stress rates are comparable to the tidal stress 
rate. The absence of a nonlinear tidal admitt- 
ance is not such a useful result for it merely 
indicates that one or more of the above con- 
ditions has not been met. 

The most interesting condition concerns the 
need for comparable tectonic and tidal stress 
rates . This condition does not arise for con- 
stitutive models of the type proposed by Stuart 
and Dieterich (1974) , nor is a nonlinear tidal 
admittance predicted. Unfortunately, tidal 
stress rates are not optimal. Mean tectonic 
stress rates are of the order 1-20 bars/year 
(10~^ - 2 MPa/year) , if it is assumed that earth- 
quake stress drops are a measure of the tectonic 
stress acciimulation between repeat earthquakes 
in the same area. The assvunptions involved in 
this estimate are discussed by Dieterich (1978) . 
Mean tidal stress rates for mid-latitudes are' 
larger, from 50-100 bars/year (5-10 MPa/yecir) , 
where the mean tidal stress rate is taken to mean 
jd^l . Consequently, nonlinear interactions are 

predicted only under what might be termed an 
accelerated tectonic stress rate. Earthquake 
precursors suggest that accelerated tectonic 
stress rates are probable before eeurthquakes . 
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Fig. 8. Possible variations in tidal admittance 
during an earthquake cycle (o^ v t) for three 

postulated earthquake models. 1 is the dilatancy- 
diffusion model in which post-earthquake dilatancy 
recovery is controlled by the diffusion of fluids. 
2 is a model in which dilatancy is recovered 
before the earthquake occurs. 3 is a model that 
postulates instantaneous dilatancy recovery at 
the time of the earthquake. The details of the 
variation of tidal admittance will depend on 
the relative tidal and tectonic stress rates, 
and the tectonic stress required for dilatancy 
to occur. 


The very appearance of a stress induced precursor 
and its. disappearance before an earthquake ' 
indicates that there is stress redistribution 
during the period for which a precursor exists. 

If precursors are a manifestation of stress 
induced dilatancy that is recovered as a broad 
zone of cracks coalesce to form a fault zone, 
there must be significant stress redistribution 
during the process. The same conclusion is 
reached if the diffusion of fluids is also 
involved. Possible variation's in tidal admittance 
for postulated pre-, co-, and post-seismic 
processes are summarized in figure 8. 

A comment on tidal triggering of earthquakes 
is relevant at this point. Heaton (1975)-, among 
others, has noted that the probability that an 
earthwake will be tidally .triggered' is greatest 
when Pg(»[d,p( . This- result- assumes that the' 

earthquake occurs through a, Mohr-Coulomb type 
failure at that phase of the tide_ when shear 
stress on the fault plane is maximized and 
normal stress is minimized. This is exactly the 
condition for which a nonlinear tidal admittance 
is least likely. Conversely, tidal triggering 
of earthquakes is least likely. when 

but this is the most favoured condition for an 
anomalous tidal admittance.' The conclusion, 
which can be tested, is that' earthquakes that 
appear to be tidally triggered are unlikely to 
have anomalous ^tidal precursors in the period 
immediately before the earthquake. .. 

The overall conclusions may be summarized in 
the following manner. The. .Earth's admittance to 
the earth tide may be sensitive to tectonic 
stress if in situ crustal rocks exhibit the same 


stress dependent properties as those observed for 
intact laboratory samples. The earth tide tilt 
and strain in a region of high and variable 
tectonic stress will exhibit linear variations if 
the Earth is nonlinear ly elastic . This result is 
similar to that predicted by Beaumont and Berger 
(1974) with the additional effect that the tidal 
admittance will be anisotropic. A more 
interesting result is predicted if crustal rocks 
exhibit stress hysteresis. At high tectonic 
stress the admittance of the Earth to the earth 
tide will be a function of the tectonic stress 
rate. In particular, a nonlinear admittance is 
predicted when the tidal and tectonic stress 
rates are approximately equal. 
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